Bistable Mott-insulator to superfluid phase transition in cavity optomechanics 
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We study the many-body state of ultracold bosons in a bistable optical lattice potential in an 
optomechanical resonator in the weak-coupling limit. New physics arises as a result of bistability and 
discontinuous jumps in the cavity field. Of particular interest is the situation where the optical cavity 
is engineered so that a single input beam can result in two radically different stable ground states 
for the intracavity gas: superfluid and Mott-insulator. Furthermore, the system we describe can be 
used as an adjustable template for investigating the coupling between cavity fields, nanomechanical 
systems operating in the quantum regime, and ultracold atomic gases. 
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Recent years have witnessed a remarkable convergence 
of interests in atomic, molecular and optical physics, con- 
densed matter physics, and nanoscience. Specific ex- 
amples include the use of ultracold atomic and molecu- 
lar systems as quantum simulators of solid-state systems 
P, the demonstration of the analog of cavity QED ef- 
fects with superconducting boxes [H, and the laser cool- 
ing of nanoscale cantilevers Q , leading to the emerging 
field of cavity optomechanics. 

The central element of most cavity optomechanical sys- 
tems consists of a Fabry-Perot type cavity with one end- 
mirror vibrating about its equilibrium position under the 
effect of radiation pressure. These devices can exhibit op- 
tical bistability, that is, the light transmitted through the 
cavity can take two distinct intensity values for a given 
incident intensity 

In this letter we show that optical bistability can lead 
to fascinating new effects in the dynamics of an ultracold 
sample of bosonic atoms trapped inside such resonators. 
In particular, at the simplest level of weak coupling and 
classical mirror motion we predict a bistable quantum 
phase transition between a Mott-insulator (MI) state and 
a superfluid (SF) state of the many-atom system. In 
the more general case where these approximations are 
removed, this system opens the way to the exploration of 
a completely new regime of interaction between light, ul- 
tracold atoms and quantum mechanical nanostructures. 

We note at the outset that clearly, a bistable transi- 
tion between a MI and a SF does not require the use of 
a cavity optomechanical system: any arrangement pro- 
ducing optical bistability would work just as well. How- 
ever, it is expected that it will soon be possible to ef- 
ficiently laser-cool one or more modes of vibration of 
moving nanoscale cantilevers or mirrors to their quan- 
tum mechanical ground state. An added advantage of 
the optomechanical cavity setup is its ability to serve as 
a diagnostic: the reflected or transmitted fraction of light 
driving the cavity has been shown to contain information 
about atomic [1, 0, [1] and mirror [i[ dynamics. This is 
what makes the coupling of ultracold atoms to optome- 
chanical systems so promising. One main purpose of this 
note is to demonstrate that these studies are rapidly be- 



coming experimentally viable. 

We recall that the MI to SF transition can occur when 
an ultracold gas of bosonic atoms is trapped by an op- 
tical lattice in the tight-binding regime Q. The ground 
state properties of the system are largely determined by 
the relative strength of the interwell tunneling energy J, 
and the intrawell pair-interaction energy U . When tun- 
neling dominates, the ground state tends to be SF. In 
the opposite case, the ground state tends to be a MI, 
characterized by a fixed atom number at each site. 

Consider then an ultracold gas of bosonic atoms 
trapped in the optical lattice provided by the standing 
optical wave inside an optical cavity exhibiting bistabil- 
ity. In the lower intensity branch the optical lattice is 
shallow, so that interwell tunneling dominates and the 
many-atom ground state is SF. In the upper branch a 
much deeper optical lattice suppresses tunneling, and the 
many-atom ground state is a MI. The state of the atomic 
system is therefore bistable, with a SF or a MI being 
formed for the same incident light field, depending on 
the history of the system. In the following we describe a 
scenario where this effect can be observed for realizable 
parameters in optomechanical resonators. 

Our study complements recent experiments on cold 
atomic gases in optical cavities with fixed ends. In each 
system two dynamical quantities are strongly coupled, 
necessitating a self-consistent, and generally nonlinear, 
description of their time evolution. Slama et al. [To| stud- 
ied the gain mechanisms behind superradiant Rayleigh 
scattering and collective atomic recoil lasing by investi- 
gating a ring cavity system. Two separate groups inves- 
tigated optomechanical systems, where collective excita- 
tions of the confined gas played the role of the mechanical 
oscillator. Brcnnecke et al. |8| demonstrated a coupling 
between a density modulated Bose-Einstein condensate 
(BEC) and the cavity field, where the phase space evo- 
lution was mapped onto that of a harmonically confined, 
cavity-coupled mechanical oscillator. Gupta et al. [hp ] 
and Murch et al. [HJ found that cavity-field coupling to 
a collective center-of-mass-motion excitation of the con- 
fined gas, resulted in oscillatory displacement of the gas. 

Our work has an especially close correspondence with 



that of Larson et al. [13j . Like us, they investigate a cold 
gas of bosonic atoms trapped by a bistable optical lattice. 
However, in contrast to our study, their cavity had fixed 
ends, and the bistability results from the strong coupling 
between the cavity field and the atomic gas. Accordingly, 
their system is modeled by a Bose-Hubbard Hamiltonian 
characterized by the parameters J, U, and chemical po- 
tential (i calculated self-consistently with the many-body 
atomic state. This self-consistent dependence results in 
a radically different ground-state phase diagram than for 
the Bose-Hubbard Hamiltonian describing our system. 
Furthermore, our system has an additional dynamical 
component - the movable end-mirror, and as already 
mentioned it can be used to investigate couplings be- 
tween these three dynamical components when one, two 
or all of them operate in the quantum regime. In this 
letter we focus on the new physics in the weak coupling 
regime and for classical mirror motion. Understanding 
this limit is an important first step in the study of the 
more complicated regimes that can be realized in our 
setup. 

On the microscopic level, the lattice potential results 
from the coupling between the intracavity field and an 
atomic resonance with frequencies, lo and u> a , respec- 
tively. As already mentioned we investigate the weak- 
coupling limit defined by Ng 2 /\A\ <C k, where N is 
the total number of atoms, go is the atom-field cou- 
pling strength, k is the cavity's natural line-width, and 
A = lu — uj a is the atom-field detuning [14| . In this limit 
the intracavity field has no significant dependence on the 
intracavity atomic population. We thus explain the gen- 
eration of the intracavity optical lattice potential by us- 
ing the theory of an empty cavity [i"o| . 

We briefly derive the necessary results from the one- 
dimensional equilibrium theory of the Fabry-Perot cavity 
shown in figure 1. The cavity consists of two mirrors, one 
fixed along x — and the other harmonically confined 
about x = Lq. Each mirror has complex transmission 
and reflection coefficients t and r, where \r\ 2 = 0.99 and 
|£|2_|_| r |2 _ ^ We only consider internal reflections where 
one may assume a 7r-phase shift, and thus we replace the 
complex r defined above with — r, where the new r is 
positive and real. The phase of t has no bearing on our 
results. 

A driving laser field E- ln of frequency lu is incident on, 
and directed normal to, the outer surface of the fixed mir- 
ror. This configuration allows a one-dimensional treat- 
ment. For a fixed cavity length, L, we follow the dis- 
cussion of Loudon [l6| to determine the transmitted in- 
tensity, /trans, exiting the cavity. The right-moving in- 
tracavity field at a; = 0, En, is determined by solv- 
ing £Jr = tE- m — rE-L under the equilibrium condition 
£l = — r exp[i2kL] En, where E^ is the left-moving cav- 
ity field at x = and k — lu/c is the wavenumber 
of the light [2(J. The resulting transmitted intensity, 
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Figure 1: Fabry-Perot cavity of length L with left-end mir- 
ror fixed along x = and right-end mirror oscillating about 
x — Lo, where L — Lo + £. The input- and transmitted- 
light intensities are labeled /in and /trans, respectively. The 
intracavity intensity at resonance (Z; n = /trans) is represented 
schematically by the sine-squared wave drawn inside the cav- 
ity. In this letter, / ~ 1 mm long, enclosing ~ 2000 standing 
wave periods, where the magnitude of the intracavity inten- 
sity iS ~ 400 /trans- 
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where E tI3lISS = tE K , I m = \E in \ 2 , and T = nr/ (l - r 2 ) 
is the cavity finesse. 

Small mirror displacements due to the intracavity ra- 
diation pressure are given by £ = ?y/trans, where rj = 
^i^ f, A is the cross-sectional area of the input laser 
beam, M is the mass of the moveable mirror and CI its 
oscillation frequency. Substituting L = Lq + £ into equa- 
tion ([T]) results in a nonlinear equation for /trans which 
is multistable with respect to /i n . We concentrate on the 
physics near cavity resonances, where kL = nir, with n a 
positive integer. For small displacements from resonance, 
the governing equation is approximately cubic in /trans, 
and predicts radiation pressure bistability [H, [poT | . 

It follows that the intracavity field intensity, I cav (x) — 
\Er(x) + El(x)\ 2 , is also bistable, and leads to a bistable 
optical lattice potential for the atoms (see Figure 2), 

VoUx) = V osc sm 2 [k(L-x)}+V h , (2) 

where P osc = ^a/ tr ans, V L = -[^ya^trans, a = 
(37rc 2 r) / (2w^A), and T is the natural linewidth of the 
atomic resonance. The microscopic origin of the propor- 
tionality constant a is the AC-Stark shift of the single- 
atom ground state. (We ignore Pl in the following since 
it is tiny compared to all relevant energies.) The position 
of the individual lattice wells is bistable as well, since a 
mirror displacement, £, displaces each optical lattice well 
by £ in the same direction. However, we consider a regime 
where £/ (7r/fe) ~ 10 -3 , and thus we ignore this effect. 

We consider a gas of ultracold bosonic atoms trapped 
in the one-dimensional optical lattice potential Vol(x). 
In the weak-coupling limit the atomic state does not alter 
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Figure 2: (Color Online) Bistability of the intracavity optical 
lattice depth, V OS c, with respect to the input light intensity 
Jin- The bistability curve is drawn with respect to the unitless 
quantities aIi n /E le and V OS c/E IC , where a = Vol (a;) /Icsvix) 
and -Ere = (h 2 k 2 ) / (2m) is the recoil energy. The detuning 
from resonance is quantified by <f>o ~ — 0.005 n, and kr\ — 
0.1 7i" a/E le , where 4>o = mod7r[feLo]. The curves AB and 
DC are the lower and upper branches of Vosc in the bistable 
region. The dashed green line connecting D and B marks the 
unstable lattice depths. The dashed gray lines DA and BC 
mark discontinuous jumps in the lattice height. 



the cavity field. Thus the atomic state is described by 
the Hamiltonian 

H= J dx^(x)(-^--^ + Vol(x) + f nO^s) , (3) 

where ip^ (x) , ip(x) are bosonic field operators, h(x) is the 
corresponding number operator, m is the atomic mass 
and g is the two-body interaction. 

We are interested in the SF-MI transition, where the 
many-atom system is accurately described by a tight- 
binding approximation that results in a single-band Bose- 
Hubbard Hamiltonian 

Hbh = - J a i a j + -n ^ ~ X ) ~ V H " 4 ' ( 4 ) 

where fit (fij) is the bosonic creation (annihilation) op- 
erator for site i, hi = aja,, and the subscript (i.j) de- 
notes a sum over nearest neighbor hopping moves. The 
tunneling matrix element is J, U is the pair interaction 
energy, and ji is the chemical potential. The param- 
eters J and U are calculated by expanding the boson 
field operators in a basis of lowest band Wannier states, 
VK 2 ') = ^-^iw{x — Xi), and then evaluating the perti- 



nent integrals 

The ground state of the many-body system described 
by equation Q is largely determined by the value of J/U, 
which depends on the intensity, wavelength and detuning 
from atomic resonance of the intracavity standing wave 
field. Of particular interest to us is the bistable regime of 
the optical lattice potential where the many-body ground 
state corresponding to the lower branch of the potential is 
a SF, while the ground state corresponding to the upper 
branch is a MI. 



Figure 3 summarizes key features of the system for a 
cavity length L — 1 mm and a moving end mirror of mass 
M = 10 mg and oscillation frequency fi = 2ir x (25 Hz). 
The cavity is loaded with a Bose-Einstein condensate of 
about 1000 sodium-23 atoms. We use an input laser of 
wavelength A = 985 nm to generate the intracavity opti- 
cal lattice potential. The optical lattice consists of about 
2000 sites, however we neglect the effects of direct atom- 
mirror interactions by assuming that only ~ 1000 sites 
near the center of the cavity are appreciably populated. 
This situation can be realized by using a gentle additional 
confining potential, and results in a lattice system with 
an average single site population near unity. 

The solid black curve in Fig. 3 is the mean- 
field SF-MI phase boundary on the axes fi/U and 



log 10 (2 J/U) [17|, [21|. For J/U smaller than this bound- 



ary the ground state is the single-particle MI, otherwise 
the ground state is SF. This diagram is overlayed by a 
plot of log 10 (2 J/U) versus I lrL for the many-body sys- 
tem described above, the logarithmic scale reflecting the 
exponential dependence of tunneling on intensity. The 
lower and upper branches are labeled with points {A,B} 
and {D,C}, respectively. These labeled points correspond 
to the lattice depths labeled in Fig. 2. For J m just above 
zero, the lattice potential is too shallow for the system to 
be described by a single-band tight-binding limit. Thus 
assigning a value of J/U is meaningless there. How- 
ever, by adiabatically increasing the intensity of the input 
laser, the condensate settles into a single band of the op- 
tical lattice potential. For high enough lattice intensity, 
the system enters the tight-binding limit. However, at 
the low-intensity edge of the bistable region, labeled A 
in Fig. 3, the system is very near the single-band tight- 
binding regime. That is, at point A a treatment with 
Eq. ^ is appropriate for our present purpose, but a 
future in-depth calculation will require including higher 
band effects. At point B, in contrast, the system is safely 
in the single-band limit, and is accurately described by 
Eq. It should be noted that the semi- log bistabil- 
ity plot has no direct correspondence to n in Fig. 3. We 
merely specify that the system is prepared so that the up- 
per branch of the bistability region lies inside the Mott- 
lobe, while the lower branch corresponds to a SF ground 
state. 

For input intensities between the points A, D (/; n ~ 
0.86 mW for our choice of parameters) and B,C (I- m ~ 
1.62 mW), the system is bistable. During an initial adi- 
abatic intensity increase, the system first resides in the 
lower branch, where the ground state is SF. Above point 
B there is only a single stable state, in the upper branch of 
intracavity intensity and lower branch of interwell tunnel- 
ing. At that point the atoms experience a much stronger 
lattice confinement, with a discontinuous phase transi- 
tion to a MI. 

The time scale over which this transition occurs is de- 
termined by the longest of the interwell tunneling times 
t ~ h/ J and the switching time of the intracavity field. 
In most cases, the intracavity field reaches a new steady- 
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Figure 3: (Color Online) Bistability of the many-body ground 
state. Thick black line: mean-field SF-MI phase boundary 
with respect to fi/U and log 10 (2J/U)). The ground state 
is the single-particle MI inside the lobe and a SF outside. 
The phase plot is overlayed with the intersite tunneling bista- 
bility curve, log 10 (2J/U) vs. Jj n curve. The (DC) branch 
corresponds to MI ground states, while the (AB) branch cor- 
responds to SF ground states. The labeled values here cor- 
respond to the lattice depths labeled in Fig. 2. The dashed 
green line indicates unstable solutions, and the arrows BC 
and DA indicate the discontinuous jumps between different 
branches. 

state value following an abrupt change in the incident 
field after a time of the order of the inverse cavity decay 
rate n ~ c\t\ 2 /L [18j]. However, for intensities switched 
from below point B to a value just above it the system 
undergoes a critical slowing down [jjjj], with a large de- 
lay before the field switches from the lower to the upper 
branch. The resulting possibility to vary the switching 
time of the light field compared to the tunneling time 



provides an important tool to investigate a variety of dy- 
namical phenomena. Most optical lattice experiments are 
performed using adiabatic tuning of the lattice height in 
an attempt to keep the system in its ground state. Af- 
ter sweeping through the discontinuity, though, we ex- 
pect that in general the many-body state will be excited 
above the ground state corresponding to the optical lat- 
tice potential. The nature of this excited state and its 
relaxation pathways are a subject of current research. Al- 
ternatively, applying a time-dependent incident field such 
that the system oscillates about the discontinuity pro- 
vides an additional tool to probe non-equilibrium prop- 
erties, and perhaps induce coupling between MI and SF 
ground states. 

Similar considerations hold when initially preparing a 
stable state in the upper branch, and then decreasing I m 
past point D, the optical lattice magnitude discontinu- 
ously jumps to its lower branch value at point A. 

In general, the setup that we described can be used to 
investigate the dynamics of coupled cold atomic gases, 
cavity fields and nanomechanical dynamics. We consid- 
ered explicitly the weak-coupling limit where the cou- 
pling between the cavity-field and the movable mirror re- 
sults in a bistable optical lattice potential for the atoms. 
We have discussed how such a cavity plus cold-atom sys- 
tem can be engineered so that SF and MI phases are 
bistable ground states for the cold-atom gas. Future 
work will extend these considerations to the situation 
where the mirror motion is quantized, and discuss in de- 
tail the dynamics of the coupled system of light, ultracold 
atoms and quantized nanostructure both in the weak and 
the strong-coupling regime. With these considerations in 
mind an important first step is to construct an experi- 
mentally viable template, where the basic physics of each 
constituent system (cold gas, cavity-field, moving mirror) 
is well understood, and one can tune the couplings. The 
setup presented above is ideal for this purpose. 
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